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Abstract. This paper reconsiders the claimed rapidity of a scheme for the 
purification of the quantum state of a qubit, proposed recently in Jacobs 2003 Phys. 
Rev. A 67 030301(R). The qubit starts in a completely mixed state, and information 
is obtained by a continuous measurement. Jacobs' rapid purification protocol uses 
Hamiltonian feedback control to maximise the average purity of the qubit for a given 
time, with a factor of two increase in the purification rate over the no-feedback protocol. 
However, by re-examining the latter approach, we show that it mininises the average 
time taken for a qubit to reach a given purity. In fact, the average time taken for the 
no-feedback protocol beats that for Jacobs' protocol by a factor of two. We discuss how 
this is compatible with Jacobs' result, and the usefulness of the different approaches. 
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1. Introduction 

Quantum feedback for continuously monitored systems has been shown to be useful for 
controlling the state of a quantum system PQ |2l El IH 13 IH] , and for improved parameter 
estimation by adaptive measurement [3 |H] . Both applications have been demonstrated 
in some prominent experiments [HI EH El Ej- Recently, Jacobs fSlE] and co-workers 
[T3] have proposed an application for quantum feedback control that incorporates 
aspects of both applications: rapid purification. As for adaptive measurements, rapid 
purification uses feedback to increase the amount of information extracted by the 
measurement. However, this information is about the final state of the quantum system, 
not a pre-existing parameter. Also unlike adaptive measurements (but as in state-control 
by feedback), the feedback in rapid purification does alter the system state on average. 
However, unlike in state-control, it is not the goal to produce a particular final state, 
only a final state of high purity (although further Hamiltonian feedback can always turn 
the final almost-pure state into any desired almost-pure state). 

In the cases considered [131 CH HHI, the objective of the feedback protocol was to 
maximise the rate of increase of the average purity of the quantum state. That is, 
to maximise the average purity in a given time. In the simplest case, a single qubit 
[T^ E], this is achieved by continually rotating the Bloch vector onto the plane that 
is orthogonal to the measurement axis. Jacobs showed that this gives rise to a factor 
of two improvement in the rate of purification over the no-feedback case. (Here Jacobs 
assumes that the measurement basis coincides with the qubit energy basis, so that in 
the no-feedback case the Bloch vector remains on the measurement axis.) 

In this paper, we demonstrate that Jacob's approach, although of undoubted 
interest and possible value, misses some important properties of the purification process. 
In particular, we show that the opposite approach — keeping the Bloch vector on the 
measurement axis — minimises the average time for the qubit to reach a given purity. In 
fact, we also find a factor of two improvement in the average time over Jacobs' protocol. 
While our results pertain only to qubits, the techniques we employ and the different 
ways to characterize purification rates in stochastic systems are potentially significant 
for broad aspects of quantum information science. 

This paper is organised as follows. In Sec. II we analyse continuous measurement 
of a qubit with no feedback, for which we use a c subscript because it has a classical 
analogue in continuous measurement of a bit. We obtain expressions for Tc, the average 
time to reach purity 1 — e, and Tc, the time for the average purity to reach 1 — e. (We 
discuss the operational meaning of the average purity in an Appendix.) In Sec. HI we 
summarize Jacobs' quantum feedback protocol for rapid purification, for which we use a 
q subscript because it has no classical analogue. We show that, for e -C 1, Tg = 2Tc even 
though Tq = Tq = Tc/2. In Sec. IV we explain how these apparently paradoxical results 
can occur by considering the (analytically derived) distributions for qubit purities under 
the different protocols. We also find numerically the full distribution of times taken to 
reach a given purity. Finally, we give a simple analytical argument for why = Tq/2, 
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which also estabhshes that the no-feedback protocol is optimal for minimizing T. We 
conclude in Sec. V with a discussion of the implications of our work for initializing 
quantum registers. We emphasize that the remarkable result that feedback allows 
an 0{d)-fo\d increase in purification rate for a d-level system, seems unlikely to be 
affected by our argument here. 



2. The No-Feedback Protocol 



2.1. Monitoring a Qubit 

We start by considering the action of continuous measurement (which can be thought 
of as continual weak measurements) on a simple model qubit, with no Hamiltonian 
(neither controlled by feedback nor constant). The conditional evolution equation for 
the qubit state matrix under continuous weak measurement of is jTHl CH Cni 

dp = dtV[a^]p + dWn[a^]p. (1) 

Here, for arbitrary operators c and p, 

V[S\p = cpc^ -^^{c'^c.p} (2) 

n[S\p = {dp - Tr[cp]p) + H.C., (3) 

while dW is an infinitesimal Wiener increment satisfying dW"^ = dt ^Tj which is the 
innovation in the measurement result jH E] . (Note that Jacobs' measurement strength 
parameter k has the value 1/2 in our units.) 

Defining z = TtI&zP] etc. we can represent p by the Bloch vector {x,y,z). By 
symmetry we can, without loss of generality, take y = 0. Then the conditional evolution 
of the qubit is given by 

dz = 2(1 - z'^)dW, (4) 
dx= - 2xdt - 2zxdW. (5) 

In terms of these variables, the purity of the qubit is given by 

p = Tr [p2] = (1 + ^ ^2^1/2. (6) 

If we choose the initial condition z = zq and x = at t = 0, then it is easy to see that 
x{t) = for all times. Then the qubit is characterized just by z{t). 



2.2. Average Purity 

Rather than using Eq. (jH), it is convenient to study the probability distribution for z 
using linear quantum trajectory theory ^3^21 • By this method, Jacobs has shown ^3] 
that for a qubit initially in the maximally mixed state {zo = 0), 2; is a random variable 
given by 



z{t) = tanh(2g). 



(7) 
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where g is a random variable having the probabihty distribution 



From this, the average qubit purity at time t is |18| ITij 




(8) 



dq. 



(9) 



cosh(2g) 



Note that we use p for the average purity, not for 1 minus the purity as Jacobs does 
[T^ . The operational significance of p is discussed in the Appendix. 

Say we are interested in reaching a certain target for the mean purity, p = 1 — e for 
e -C 1. This will occur in the long time limit t 3> 1, but even for t = 5, it is a very good 
approximation to replace Eq. with the asymptotic {f^V) expression 



J-oo cosh(2g) 4^27^ 
Thus, the time t = at which the average purity reaches the level 1 — e is, in the 
asymptotic regime of e <^ 1, 



Here the subscript c denotes classical, because this measurement involves no coherences 
(it applies equally to a measurement of a classical bit as to a qubit). This is because, 
as noted above, the action of the measurement of Oz does not move the Bloch vector 
away from the z-axis and we are ignoring any Hamiltonian evolution. In more general 
situations, where the intrinsic qubit Hamiltonian does not commute with o"^, the Bloch 
vector will move away from the measurement axis. If this is the case, Hamiltonian 
feedback control would be required to rotate the Bloch vector back onto the measurement 
axis in order to reduce the problem to the classical one solved here. 

2.3. Average Purification Time 

Rather than the time taken for the average purity to reach a given level, we may be more 
interested in the time at which a particular system reaches a target purity p = 1—e. For 
this it is more convenient to return to the (nonlinear) quantum trajectory equation (j3)). 
From this equation, an ensemble of qubits with z = zq and x = at t = is represented 
by the probability distribution p{z;t\zo;0) obeying the Fokker-Planck equation jT7] 



with initial condition p{z] 0) = 6{z — Zq). We wish to study the time T at which \z\ first 
attains the value Z such that (1 + Z^) /2 = 1 — e. We can do this by considering the above 
Fokker-Planck equation with absorbing boundary conditions at Z and —Z. Writing its 
solution as p{z, t) [where the tilde denotes the presence of boundary conditions], we can 
thereby define the probability that T > t for an initial value of zq as ^\ 




(10) 






(12) 




(13) 
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It is simple to see that the mean time until the qubit is purified to the desired level 
is given by 

T{zo) = / G{t\zo)dt (14) 
Jo 

Following the method in Ref. this quantity (known as the average time of first 
passage) obeys the ordinary differential equation 

2(1 - z',Yr'izo) = -1, (15) 

where the dash represents differentiation with respect to zq- The boundary conditions 
are T(— Z) = T{Z) = 0, since a qubit with \zo\ = Z is at the boundary already so T = 
there by definition. The solution to this problem is 

T{zq) = -(ZartanhZ — ZoSntanhzQ) . (16) 

We are again interested in the case where the qubit is initially in a completely 
mixed state: Zq = 0. Then T = (l/4)ZartanhZ. For e small we have Z ~ 1 — e, and 

f,^iln(2/e)~iln(e-i) (17) 

(Note the subscript c again.) That is, the mean time Tc to reach a purity of 1 — e is one 
quarter the size of the time to reach a mean purity of 1 — e. 

3. Jacobs' Protocol 

Jacobs' protocol is to maximize the increase in the average purity at every point in time. 
From Eqs. (j^l)-® it is easy to show (using the Ito calculus [Ej) that 

ds = -{8s^ + Ax^s)dt - AzsdW. (18) 

Here s = 1 — p is sometimes known as the linear entropy. Note the distinction between 
s, a time-dependent random variable, and e, a fixed parameter related to a target purity 
to be achieved by a protocol. It is clear that on average s will decrease most rapidly if 
X is maximized. That is, if one uses feedback to rotate the Bloch vector onto the x axis. 
Then one has 2 = and x"^ = 1 — 2s, so 

ds = -Asdt (19) 

Thus, under Jacobs' rapid purification scheme the evolution under ideal Hamiltonian 
feedback control is deterministic! The purity of the qubit is given by 

Tr[p^] = 1 - \e~'\ (20) 

identically for all qubits, so p = p, and also T = T = r. That is, the time of first passage 
and the time for the average purity to reach the desired level are equal. For a target 
purity of 1 — e, these times are 

1 

4 4 

Here the subscript q stands for quantum, since this adaptive technique exploits the qubit 
coherences and cannot be applied to a classical bit. For e <^ 1, is exactly twice as 
long as Tc, even though Tg is half as long as Tc. 



= Tg = --ln2e~ -ln(e"^) (21) 
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4. How is this possible? 
4.I. Distributions of Purity 

Some insight into these results can be found from considering p{p; t)dp, the probability 
distribution for the purity at time t. In terms of this, the mean purity at a given time 
t is p{t) = Jo pp{p]t)dp, so the time r at which the purity reaches the desired value is 
defined as 

l-e= !\p{p;T)dp (22) 





To define T it is necessary to consider a different distribution, p(p; T)dp, which arises 
from the absorbing boundary conditions as described in Sec. II. In terms of this, the time 
T is given by 

_ /"OO rl—e 

T = dt dpp{p;t). (23) 



10 JO 

(This is derived in the same way as Eq. (|Tl|l .) However, when e is small it is very unlikely 
(with probability of order e in fact) for the purity of the qubit to reduce significantly 
below 1 — e once it has crossed that barrier. Thus p{p; t) is sufficiently similar to p{p] t) 
that it can help us to understand T also. From Eqs. (jHl)^®, it follows that for the 
classical (no-feedback) measurement the distribution of purities is 

p-2t exp ( — ■5Tartanh^-v/2iD — 1 ) 
p^ip; t)dp = —= '^-^==Z=^J-dp. (24) 

By contrast, under Jacobs' quantum feedback protocol 

p,{p-t)dp = 5{p-l + e-'^'/2)dp (25) 

The two purity distributions are plotted in Fig. 1 for two different times, t = Tc and 
t = Tq = Tq. With e = 10"^ these times are t = f^ = 1.72694 and t = Tq = = 3.2806. 
The agreement between the numerical histograms and theoretical results are good; the 
slight differences between the two are related to the limited numbers of runs {circa 20 
000) and rounding errors in the numerical calculations for very low e values. In each 
case. Fig. 1(a) and 1(b), the majority of the probability distribution for the classical 
(stochastic) scheme lies to the right of the quantum (deterministic) scheme which is 
indicated by the red line. This means that the majority of cases should produce a 
significantly lower value for the purity than Jacobs' approach, whilst the average purity 
for the stochastic case is unduly affected by the relatively few instances where the purity 
is still relatively far from one. 



4.2. Distributions of Purification Time 

In Fig. 2 we plot the classical and quantum probability distribution for the time of first 
passage (T), as determined from numerical simulation of (0)) and (for Jacobs' protocol) 
from Eq. ()21|) . for the same value of e. As we would expect from Fig. 1, the majority 
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Figure 1. Probability distribution of purity values generated by the no- feedback 
scheme — theoretical (black dashed line) and numerical (blue histogram) — and 
Jacobs' deterministic scheme (red solid line) at different times: (a) t = 1.72694 and 
(b) t — 3.2806. Note the logarithmic scale for s = 1 — p. 



'(/} 



CL 




Time of First Passage 



Figure 2. Probability distributions of average times of first passage for e = 10 ^ for 
present scheme (blue histogram) and Jacobs' deterministic scheme (red solid line). 
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Figure 3. (a) Average time to first passage for Jacobs' sclieme (solid red line) and the 
scheme presented here (solid blue line) with dashed lines indicating standard deviation 
of the first passage distribution for Tc, (b) Ratio of average first passage times Tq/Tc- 
data joined by solid line; theoretical curve ln(l/2e)/Vl — 2eartanhVl — 2e (black 
dashed line which asymptotes to 2 and is partially covered by numerical results). Both 
graphs are plotted as functions of e on a log scale. 

of classical (stochastic) cases reach the desired purity level before they would using the 
deterministic, ideal feedback scheme. In Fig. 3(a) we plot the average first passage time 
as a function of e, and in Fig. 3(b) the ratio of average first passage times produced by 
two measurement schemes as a function of e. Fig. 3(b) shows that the relative benefit 
of the classical approach tends to the theoretical value of Tq/T^ — 2 as e — 0. Note 
that the improvement offered by our classical protocol is significant even for quite low 
target purities (e.g. e = 0.01). Moreover, Jacobs' purification time and the average 
first passage time are separated by at least one standard deviation of the first passage 
time for all e < 0.003. These facts are encouraging for experimentally demonstrating 
the differences between the schemes, even if Jacobs' ideal Hamiltonian feedback scheme 
cannot be realized perfectly in the laboratory. 

In Jacobs' scheme, the purity of the state always increases. The observer has no 
prior knowledge about the observable dz being measured because the state is always 
rotated to lie on the horizontal plane of the sphere. As a consequence, any result always 
increases the observer's certainty about the state. 
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By contrast, in the classical scheme analysed here, the measurement may increase 
or decrease the observer's certainty about the system. Consider the situation where the 
qubit is in a mixture of being in state o"2 = 1 and o"z = —1, but more likely to be in the 
former so that ^ > 0. From Eq. if dW > in the next infinitesimal interval, then 
the observer will become more certain that (Tz = 1 is the true state of the qubit, and z 
(the observer's expectation value for a^) will increase. On the other hand, if dW < 0, 
the observer will become less certain and z will decrease. That is, in any infinitesimal 
interval, the purity decreases half the time! Why then does the purity increase at all? 
The answer is that p varies as z"^, so even though the positive and negative changes in z 
balance, the change in p is positive on average. But, not surprisingly, the increase in p 
is slower than under Jacobs' scheme. Also, in the classical case there is clearly going to 
be a far greater spread in the purity than in Jacobs' scheme, as shown in Fig 1. But this 
is exactly what allows for the average time T to reach a purity threshold to be shorter 
in the former case. 

4-3. Simple Derivation and Optimality 

The results already presented in this section confirm that the no-feedback protocol gives 
an improvement over Jacobs' protocol for the average time to purify a qubit, and give 
some insight into how this is possible. However, what is lacking is (a) a simple derivation 
of the factor of two improvement, and (b) a demonstration that the no-feedback protocol 
is optimal for this purpose. We now provide both of these. 

By considering the average of Eq. (fTHj) it was easy to see that Jacobs' protocol 
maximized the decrease in the average linear entropy s = 1 — p. However, as argued 
above, for the no-feedback case this average gives undue weight to the tail of the 
distribution in which s is much larger than the mean. This can be rectified by considering 
Ins, the distribution of which is much more symmetrical, as seen in Fig. 1. From 
Eq. (fTH|) . the stochastic differential equation for Ins is, using the Ito calculus [T7], 

d\ns = -4:{2s + x^ + 2z^)dt-4:zdW (26) 

Taking the average of this equation, we see that the most rapid purification (for a given 
s) is when z"^ is maximized at 1 — 2s and set to zero; that is, exactly the classical 
protocol considered above. Moreover, in this classical case we obtain, for almost-pure 
states (z^ ~ 1), 

rf(lns)/dt ~ -8, (27) 

which is the characteristic rate for purification to a high level of purity. By contrast, 
under Jacob's protocol we set z"^ = and = 1 — 2s to obtain 

d\ns/dt = -4:. (28) 

This is of course the rate of purification found in Eq. (fTn|l . and is slower than the classical 
rate ^T7\i by the factor of two found previously. 



Reconsidering Rapid Qubit Purification by Feedback 



10 



5. Conclusions 

In this paper, we have considered some interesting properties of the rapid purification 
protocol proposed by Jacobs for quantum systems subject to continuous measurements 
and Hamiltonian feedback control. In particular, we have demonstrated that, although 
Jacobs' scheme (keeping the Bloch vector perpendicular to the measurement axis) 
provides the fastest increase in the average purity of the qubit, the opposite approach 
(keeping the Bloch vector on the measurement axis) provides the shortest average time 
to reach a given purity level. This counter-intuitive result is due to the following: 
in Jacobs' approach the purification is deterministic, whilst the classical approach is 
stochastic where the distribution of purities at a given time is heavily skewed. Thus, 
although the majority of qubits will reach a given purity level quicker than under Jacobs' 
scheme, there is a minority of qubits that have a relatively low purity, which reduces 
the expected (average) purity below that of Jacobs' scheme. 

Aside from their theoretical interest, these results have most obvious application in 
the initialization of qubits for quantum information processing. The fact that, under the 
classical protocol, half the qubits end up in the = 1 state and half in the = —1 state 
does not matter; it is known which state the qubit is in so it can always be rotated to 
the desired fiducial state. Since it is increasingly common to consider non-deterministic 
protocols in quantum computing, it seems likely that a scheme that minimizes the 
average time to purify the qubits in an ensemble would be preferred. The fact that no 
Hamiltonian feedback control is required for the classical protocol is another advantage 
it has over Jacob's protocol. However, if it were necessary to purify all members of an 
ensemble of qubits to a given degree of purity in as short a time as possible, then Jacobs' 
rapid purification by quantum feedback should be used. 

Finally, we wish to emphasize the following. In very recent work. Combes and 
Jacobs ^H] have shown that Hamiltonian feedback control can be used to purify a qudit 
(that is, a d-level quantum system) a factor of order d faster than can be achieved 
without feedback. There is nothing in our analysis to suggest that this factor of d would 
change if one calculated the average time to purify (as we considered here) rather than 
the average purity, as Jacobs has always considered [131 CH HSj . We suggest that any 
difference between the two types of purification time would be of order unity, but that 
remains to be confirmed by future work. 
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Appendix 

One could question wlietlier tlie average purity p as used by Jacobs is a meaningful 
concept operationally. If one were to monitor a qubit (or bit) as above, then forget 
the particular measurement record obtained, and then ask what is the purity of the 
resulting state, the answer is not given by Eq. 0. The reason is that purity is a 
nonlinear function of the state, unlike z for example. Of course if one were to ignore 
completely the measurement record then the purity of the resultant state would be 1/2, 
because the state is as likely to be purified towards z = lastoz = — 1. Thus it is only 
sensible to say that at the end of the measurement the state is first rotated (if necessary) 
so that 2; > 0, and then the measurement record is forgotten. The resultant state is 
then defined by the expected value of cr^, which in terms of the unrotated stochastic 
variable z is given by From above, this evaluates to 

(1^1) = 2 / smh{2q)p{q)dq = 1 - / e'"^ l^'-'^Uq. (29) 

h \/27Tt Jo 



/27rt 

The purity of this average state is Pav. = 1 + (kl)^- Putting = 1 — e, in the limit 
t 1 this evaluates to 

' e-^''dq=^=. (30) 



V27rt ^0 V27rt 
This is only a constant factor (4/7?) larger than that calculated above (|l(Jj). so the 
asymptotic expression for in Eq. is unchanged. For these reasons, in the body of 
this paper we just use Jacobs' p for simplicity. 
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